Measurement of the Electric Form Factor of the Neutron at Q"^ = 0.5 and 1.0 GeV^/c^ 
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The electric form factor of the neutron was determined from measurements of the d{e,e'n)p 
reaction for quasielastic kinematics. Polarized electrons were scattered off a polarized deuterated 
ammonia (^^NDs) target in which the deuteron polarization was perpendicular to the momentum 
transfer. The scattered electrons were detected in a magnetic spectrometer in coincidence with 
neutrons in a large solid angle detector. We find ~ 0.0526 ± 0.0033{stat) ± 0.0026(sj/s) and 
0.0454 ± 0.0054 ± 0.0037 at = 0.5 and 1.0 (GeV/c)^ respectively. 

PACS numbers: 14.20.Dh, 13.40.Gp, 24.70.+S, 25.40.Ve 



The electric form factor of the neutron is a fun- 
damental quantity in nuclear physics. Knowledge of 
over a large range of momentum transfer leads to insight 
to physics beyond the simple SU(6) symmetric models, 
for which it must vanish. Beyond nuclcon structure, our 
understanding of G^ has an impact on a broad range of 
topics, which vary from the charge radii of nuclei 1] to 
extracting the strangeness content of the nucleon . 

Historically, measurements of G^ have suffered from 
large uncertainties due to experimental technique and nu- 
clear model dependence. Early attempts to measure G^ 



from cross sections sometimes failed to determine even 
the sign. Until recently, the best determination of G^ 
came from elastic electron-deuteron measurements, but 
the errors were large, ^30%, due to their dependence on 
the nuclcon-nucleon potential model j^J. The first po- 
larization measurements of G^, conducted with ^H and 
■^He targets, differed significantlv because final state in- 
teractions were not addressed 0, 0| . These experiences 
highlight the importance of measuring G^ using different 
reactions. 

In the last few years the experimental understanding 
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of has improved considerably. The results from po- 
larization experiments described in Ref. are con- 
sistent and provide good accuracy in the kinematic re- 
gion of the four-momentum transferred squared < 
0.7 (GeV/c)^ (henceforth units of are assumed to be 
(GeV/c)^). Until now, values of determined from the 
deuteron quadrupole form factor [l^ provide the only 
information in the kinematic region above > 0.7 and 
leave the > 1.6 region undetermined. 

This letter describes a first measurement of G^ at 
= 1.0 using a polarized target. In addition, the result 
for ~ 0.5 has half the uncertainty as the previous mea- 
surement (3] and is the most precise result near the peak 
of G^. These results are largely insensitive to the model 
of the nucleon-nucleon (NN) potential so that, compared 
to those of Ref. they are more reliable. 

To determine G^, the helicity dependent rate asymme- 
try in electron scattering was measured. In the ideal case 
of a polarized electron scattering elastically off a free po- 
larized neutron, with the neutron polarization vector in 
the scattering plane and perpendicular to the momentum 
transfer q, Gp is related to the beam-target asymmetry 
term A^^ ilj by 

,v ^ -2VT(r + l)tan(0j2)G^G;} 

(Gs")2+T[l + 2(l + T)tan2(0,/2)](Gi})2' 

where r = Q'^/AM^, Mn is the mass of the neutron, G/} 
is the magnetic form factor of the neutron, and de is the 
electron scattering angle. 

For lack of a free neutron target, the actual measure- 
ments were performed on a polarized deuterium target. 
The measured experimental asymmetries were due to a 
combination of several physics asymmetries scaled by the 
electron and target vector and tensor polarizations (that 
are described in detail in our previous work 0, not re- 
peated here for brevity's sake). The proper averaging of 
the asymmetry (symmetrically around q) and the negligi- 
ble contributions from the target tensor asymmetry sim- 
plify the relationship of the measured asymmetry e to the 
deuteron vector asymmetry A^^ so that e = fPePiA^^, 
where Pe is the beam polarization, P^ is the deuteron 
vector polarization and / is the dilution factor due to 
scattering from nucleons other than polarized deuterons 
in the target. Calculations show that depends lin- 
early on G£ for the kinematics of the experiment . 

The measurements were conducted in Hall C of the 
Thomas Jefferson National Accelerator Facility in a setup 
similar to that of the previous measurement 0- The 
longitudinally polarized electron beam was scattered 
off a polarized frozen deuterated ammonia (^^NDa) tar- 
get. The scattered electrons were detected by the High 
Momentum Spectrometer (HMS), and the neutrons were 
detected by a dedicated neutron detector. The central 
kinematics for the two measurements as well as the aver- 
age beam and target polarizations are listed in Table 
The average deuteron luminosity was lO'^^/cm^/s. 
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TABLE I: Central kinematics and average polarizations. E 
{E') is the energy of the incident (scattered) electron. 6^ (On) 
is the angle of the scattered electron (neutron). 6b is the 
orientation of the target polarization axis. (Pe) and (Pt) are 
the average beam and target polarizations, respectively 

The polarized target consisted of ammonia gran- 
ules submerged in liquid He and maintained at 1 K by 
a ^He evaporation refrigerator. The target spins were 
aligned by a 5 T magnetic field generated by a pair of su- 
perconducting coils. The polarization was enhanced via 
dynamic nuclear polarization Il7l and measured with a 
continuous- wave NMR system To minimize local- 

ized heating and depolarization of the target material, 
the electron beam was rastered uniformly in a circular 
pattern with a 1 cm radius. 

A two-magnet chicane compensated for the deflection 
of the electron beam by the target field. After traversing 
the target, the beam passed through a helium bag to a 
special dump in the experimental hall. The polarization 
of the beam was measured at regular intervals through- 
out the experiment with a M0ller polarimeter 0|. The 
beam helicity was changed in a pseudorandom sequence 
30 times per second to minimize charge fiuctuations cor- 
related with helicity. 

The HMS was operated in its standard mode for the 
detection of electrons. The established reconstruction 
algorithms were augmented to account for the large beam 
rastering and the effects of the target magnetic field on 
the scattered electrons. 

The neutron detector consisted of multiple vertical 
planes of segmented plastic scintillators. Two planes of 
thin scintillators served to distinguish charged particles. 
Behind these were six planes of thick scintillators to de- 
tect the neutrons. All scintillators were equipped with 
photomultipliers on both ends to provide spatial and time 
information for the detected particle. The 88 thick scin- 
tillators provided a neutron detection volume that was 
roughly 160 cm wide, 160 cm tall and 90 cm deep. The 
front of the detection volume was approximately 4.2 m 
and 6.2 m from the target during the = 0.5 and 1.0 
measurements, respectively. The detector was shielded 
from direct gamma rays from the target by a 2.5 cm lead 
curtain, and the entire assembly was housed in a thick- 
walled concrete hut, which was open to the target. 

The trigger was set up so that the neutron detector was 
read out for every electron trigger in the HMS. Coinci- 
dences between the electrons and the knock-out nucleon 
were determined offline. 

The experiment was simulated using Monte Carlo 
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FIG. 1: Comparison of = 1.0 MC and data electron- 
neutron event distributions from all materials in the target 
for four kinematic variables: invariant mass W , Q^, angle be- 
tween the neutron and q in lab frame Onq, and angle between 
neutron-proton system and the momentum transfer in the 
center-of-momentum frame 9"^. The solid grey histograms 
correspond to the data, and the dotted black histograms cor- 
respond to the simulation. Only the shaded regions, domi- 
nated by events from were used to determine . 



(MC) techniques. The simulation package, based on 
MCEEP ^2^, included the charged particle transport 
through the target's magnetic field and the optical and 
aperture model of the HMS from the Hall C simulation 
package SIMC j2l|. The MC served two principal func- 
tions: to determine the dilution factor and to average 
the theoretical asymmetries over the acceptance. For the 
dilution factor, all target materials were included in the 
simulation: the deuterium and the nitrogen in the am- 
monia, the liquid helium bath, the aluminum target cell 
windows and the NMR coil. 

Contributions from electron-neutron events originating 
from quasielastic scattering and pion production were in- 
cluded in the MC for all target materials. It was found 
that the (e, e'rCji: contribution to the event sample was 
neghgible for = 0.5 and less than 0.5% for = 1.0. 
Two-body knockout contributions were examined in the 
MC (following jl^l) and were also negligible. 

The simulated rates were normalized to the measured 
ammonia rates, and the variations in the ratios of MC 
rate to observed rate for various target materials were 
used to determine the uncertainty in the dilution factor. 
A comparison of event distributions of the data and MC 
for several kinematic variables for = 1.0 is shown in 
Fig. n agreement is very good. The comparison for the 
= 0.5 data is similar and slightly better than in the 
previous experiment . 



Several corrections to the measured asymmetry were 
applied. The leading correction was for charge ex- 
change of protons in the lead shielding: —3.8 ± 1.1% 
(-3.3 ± 3.0%) for Q2 = o.5 (Q^ = i.o). The charge ex- 
change in the target material itself was small: — 0.3±0.3% 
(—0.3 ± 0.3%). The accidental background rate was 
1.9 ± 0.1% (0.5 ± 0.2%) with no statistically significant 
asymmetry. Internal radiation effects on the asymmetry 
were also corrected: 0.9±0.5% (0.6±0.5%). The effect of 
external radiation was insignificant compared to the 0.5% 
MC statistical uncertainty. The contamination from neu- 
tral pions generated by protons in the lead shielding of 
the neutron detector was also found to be insignificant. 

The physics model of A^^ used in the MC was based 
on the calculations of Arenhovel, Leidemann and Tomu- 
siak 01 . It included a non-relativistic description of the 
n — p system in the deuteron using the Bonn i?-space 
NN potential for both the bound state and the final 
state interactions. The full calculation included meson 
exchange, isobar configuration currents and relativistic 
corrections. 

The model assumed a scaled Galster parameteriza- 
tion for and the dipole parameterization for G^" . 
The Galster parameterization is G^{Q'^) = — ^„t/(1 + 
5.6r)Gi5(Q2), where Go = 1/(1 + QV0.71)2 is the dipole 
form factor. Various scale factors of this parameteriza- 
tion, ranging from 0.5 to 1.5, were examined. The po- 
tential impact of the dependence of the G£ param- 
eterization was found to be negligible because the 
acceptance was not large and the dependence of G£ 
over the acceptance was mostly linear. The acceptance- 
averaged value of A^^ has a linear dependence on G^}, so 
one can trivially correct for more accurate G^" values. 

The value of G^ was determined by comparing the ac- 
ceptance averaged of the data to that of the MC. The 
theoretical asymmetries were determined for a range of 
scaling factors of the Galster parameterization to deter- 
mine the corresponding G£: G^YGalster — 1.003±0.064 
and 1.172 ± 0.140 for = o.5 and 1.0, respectively. 
The uncertainties are statistical only. To account for the 
slight deviations of G^ from the dipole form factor, the 
recent fit to the G^ data [l^ was used: G^J^nGo = 
1.007 ± 0.005 and 1.072 ± 0.014 for = o.5 and 1.0, 
respectively. Then the values for G^ are: 

G^(g2 = 0.5) 0.0526 ±0.0033 ±0.0026, 
G^(g2 ^ 1.0) = 0.0454 ±0.0054 ± 0.0037, 

where the first uncertainty is statistical and the second 
is systematic. The = 0.5 result agrees well with the 
previous result reported in Ref. j^. 

Other systematic uncertainties for = 0.5 (Q^ = 
1.0) were: dilution factor, 3.4% (3.0%); target polariza- 
tion, 2.9% (4.6%); central kinematic values, 1.2% (3.4%); 
beam polarization, 1.1% (3.3%). 

The results as compared to recent measurements are 
shown in Fig. |21 The new data at = 1.0 provides the 
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FIG. 2: Comparison of this experiment with data from recent 
measurements. The data points are: diamonds (T^. hoUow 
squares 10], triangle [J, hollow circle 0, solid square 0, plus 
sign 11]. Solid circles: this paper. See text for a description 
of the curves. 

important experimental confirmation of the decline of 
following the Galster form. The data shown in the figure 
were fit to the form G£{Q'^) = — /i„aT/(l + pt)Gd{Q^), 
The parameter a = 0.895 ± 0.039, the slope at ^ q, 
is fixed by atomic measurements of the neutron charge 
radius Hg. The fit yielded p = 3.69 ± 0.40. The errors 
obtained for the fit parameters are uncorrelated. The 
one-sigma error region of the fit is shown in Fig. [21 as the 
shaded band. 

Many recent models [23|-|13| have attempted to 
predict or fit the nucleon electromagnetic form fac- 
tors. Fig. 121 compares the data with two recent calcu- 
lations that use covariant formulations of the constituent 
quark model with quark-quark interactions fitted to spec- 
troscopic data. The point-form spectator approxima- 
tion (PFSA) of 113 (dot-dash) uses a Goldstone boson 
exchange interaction with pointlike constituent quarks 
while the light-front (LF) calculation of 28] (short dash) 
uses a one-gluon exchange interaction with constituent 
quark form factors fitted to data for < 1. The use of 
constituent form factors improves the fit to the nucleon 
magnetic form factors at larger , but the PFSA seems 
to describe G^ better at low with fewer parameters. 

Also shown in Fig. (21 are the results from a hy- 
brid model that interpolates between vector-meson dom- 
inance at low and perturbative QCD at high 
113 (solid line), from a light front model where the 
nucleon is considered a system of three bound quarks 
surrounded by a cloud of pionslsjj (long dash) and from 
a soliton modelji^ (dotted) whose basic features include 
an extended object, partial coupling to the vector mesons 
and relativistic recoil corrections. While all these models 
agree qualitatively with the data, none agree with the 
data for the entire range of Q^. 

It is remarkable that in the last five years, the exper- 
imental precision in G^ measurements has improved to 
the 10% level. This significant improvement provides a 



rigorous challenge for models of the nucleon structure 
because this electromagnetic form factor is the most sen- 
sitive to physics beyond the simplistic SU(6) symmetric 
picture. 

In conclusion, the electric form factor of the neutron 
at — 0.5 and 1.0 has been determined from measure- 
ments of the beam-target asymmetry. This experiment 
provides the highest measurement to date using a po- 
larized target and the most precise measurement near the 
maximum of G^ . 
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